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Abstract 

We report on the observation of B — > Z}*°7r + 7r + 7r _ 7r~ decays. The branching 
ratio is (0.30±0.07±0.06)%. Interest in this particular mode was sparked 
by Ligeti, Luke and Wise who propose it as a way to check the validity of 
factorization tests in B —> D* + ir + ir~ ir~ ir° decays. 
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Factorization is the assumption that in two-body hadronic B decays the decay amplitude 
can be expressed as a product of two currents, just as in semileptonic decays where one 
current is hadronic and the other leptonic. Use of factorization has been crucial in creating 
models for understanding the underlying weak decay dynamics 0. 

In previous work we found a large branching fraction of (1.72±0.14 ±0.24)% for the decay 
E? -> D*+ 7T + 7r 7T 7T° H]. This reaction can proceed via several possible tree level processes. 
The simplest diagram, shown in Fig. |l](a), has the four-pions emitted from the virtual W~ . 
Assuming that this is indeed the dominant process, Ligeti, Luke and Wise (LLW) [[| have 
compared the 47r~ invariant mass spectrum with r~ — > 7r + 7r~7r~7r°z/ data [[|. Using a model 
based on factorization they show that the data agree up to a 47r~ mass-squared of 2.9 GeV 2 , 
within a precision of about 15%. 

However, the agreement may be fortuitous, rather than a success of factorization, if other 
diagrams are present. For example another possible diagram is shown in Fig. 0(b), where the 
D* + and the ir° are produced at the lower vertex and the virtual W~ manifests as 7r + 7r~7r~. 
This process was searched for in the original publication. Definite evidence was lacking but 
a stringent upper limit could not be set. 

Here we search for the process B — > D* 7t + 7t + 7t~tt~ as suggested by LLW. This can 
be produced by the diagram in Fig. 0(b), where the D*° combines with one of the 7r + 's to 
form a low-mass system. It could also be produced by the color-suppressed process shown 
in Fig. |](c). In this paper we indeed show that the process B — > D* Q n + n + n~n~ has a 
significant branching ratio and try to ascertain the dominant production mechanism. The 
data sample consists of 9.0 fb _1 of integrated luminosity taken with the CLEO II and II. V 
detectors || using the CESR e + e~ storage ring on the peak of the T(4S') resonance and 
4.4 fb _1 in the continuum at 60 MeV less center-of-mass energy. The sample contains 19.4 
million B mesons. 

Hadronic events are selected by requiring a minimum of five charged tracks, total visible 
energy greater than 15% of the center-of-mass energy, and a charged track vertex consis- 
tent with the nominal interaction point. To reject noia-BB continuum we require that the 
Fox- Wolfram moment i?2 be less than 0.3 ||. Track candidates are required to pass through 
a common spatial point defined by the luminous region. Tracks with momenta below 900 
MeV/c are required to have an ionization loss in the drift chamber within 3 standard devi- 
ations of that expected for their mass hypothesis (tt/K). Photon candidates are required to 
be in the "good barrel region," within 45° of the plane perpendicular to the beam line that 
passes through the interaction point, and have an energy distribution in the Csl calorimeter 
consistent with that of an electromagnetic shower. To select 7r°'s, we require that the dipho- 
ton invariant mass be between —3.0 to +2.5<r of the 7r° mass, where the resolution a varies 
with momentum and has an average value of approximately 5.5 MeV. The 7r° candidates are 
then kinematically fit by constraining their invariant mass be equal to the nominal ir° mass. 

We select D° candidates in the K~tt + decay mode. We require that the invariant mass 
of the D° candidates lie within ±2.5<r of the known D° mass. The D° mass resolution 
varies with D° momentum, p, as px0.93xl0~ 3 +6.0 (units of MeV). We use the analogous 
requirement on the D*°-D° mass difference. In this case the mass difference resolution is 
0.90 MeV. 
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FIG. 1. Diagrams for B — ► D*Trinrir decays, (a) Charged current tree level diagram for 
D* + (47r)~ . (b) Charged current tree level diagram for (D* + tt°) 7r + 7r~7r~ or (D*°tt + ) tt + tt~it~ . The 
D*tt system can form a D** resonance, (c) Color suppressed diagram for D*°(47r)°. 

We start by looking for the D*°(4ir) final state.Q The D*° candidates are pooled with 
all combinations of 7r + 7r + 7r~7r~ mesons. Next, we calculate the difference between the beam 
energy, E beam , and the measured energy of the five particles, AE. The "beam constrained" 
invariant mass of the B candidates, M B , is computed from the formula M B = E$ eam — 
(Si Pi) 2 - To further reduce backgrounds we define 

2 = / am p . y / m(kk)-m d y ( m( 77 )-m.o y 

Xb \a(AM D *) J + \a(M(Kn) - M D ) J + ^( M (77) ~ M n o) J ' U 

where AMp. is the computed D*° — D° mass difference minus the nominal value and the 
cr's are the measurement errors. We select candidate events requiring that x% < 5. 

We show the candidate B mass distribution, M#, for AE in the side-bands from —6.0 
to —4.0a and 4.0 to 6. Oct on Fig. |2|(a). The AE resolution is 14 MeV (a). The sidebands 
give a good representation of the background in the signal region. We fit this distribution 
with a shape given as back(r) = pir^/l — r 2 e~ P2<yl ~ T \ where r = M^/5.2895, and the Pi are 
parameters given by the fit. 

We next view the Mb distribution for events having AE within 2a of zero in Fig. |2|(b). 
This distribution is fit with a Gaussian signal function of width 2.7 MeV and the background 



In this Letter (4-7r)° will always denote the specific combination tt + tt + tt tt 
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FIG. 2. The B candidate mass spectra for the final state D*°7r + 7r + 7r~7r~, (a) AE sidebands, 
(b) for AE consistent with zero. 

function found above whose normalization is allowed to vary. The Gaussian signal width is 
found from Monte Carlo simulation. The largest and dominant component results from the 
energy spread of the beam. We find a total of 64±16 events, thus establishing a signal. 

The error due to the background shape is evaluated in three ways. First of all, we change 
the background shape by varying the fitted parameters by la. This results in a change of 
±9.3%. Secondly, we allow the shape, p 2 , to vary (the normalization, pi, was already allowed 
to vary). This results in 11% decrease in the number of events. Finally, we choose a different 
background function, back'(r) = p\r\f\ — r 2 (1 + p 2 r + p 3 r 2 + p 4 r 3 ) , and repeat the fitting 
procedure. This results in a 9.3% increase in the number of events. We assign ±7 event 
systematic error due to this source. 

We have investigated two mode-specific backgrounds that could, in principle, induce fake 
signals. These include the final states D* + 7t + 7t~7t~tt where we miss the slow tt + from the 
D* + decay and the 7r° and D° happen to satisfy the D*° requirement, and .D°7r + 7r + 7r~7r~7r°, 
where the D° and the ir° happen to satisfy the D*° requirement. We find that the efficiencies 
for each of these modes to contribute are small. The first final state was measured as having 
a branching ratio of 1.72% 0. It would contribute 0.4±0.5 events. The second final state has 
never been measured. It would contribute 1.6 ± 0.5 events per 1% branching ratio. Taking 
into account all sources of systematic error we observe 64±16±7 signal events. 

In order to find the branching ratio we use the Monte Carlo-determined efficiency, shown 
in Fig. |3] as a function of (4tt)° mass. Since the efficiency varies with mass we need to deter- 
mine the (47r)° mass spectrum in order to determine the branching ratio. To rid ourselves 
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of the problem of the background shape, we fit the B candidate mass spectrum in 100 MeV 
bins of (47r)° mass. The resulting (47r)° mass spectrum is shown in Fig. f|. We find 
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FIG. 3. The efficiency for the final state Z)* 7r + 7r + 7r~7r 



B(W -»• D*\ + 7r + Tr'7r~) = (0.30 ±0.07 ±0.06)% . (2) 

The systematic error arises mainly from our lack of knowledge about the tracking and 7r° 
efficiencies. We assign errors of ±2.2% on the efficiency of each charged track, and ±5.4% 
for the 7T°. The total tracking error is found by adding the error in the charged particle track 
finding efficiency linearly for the 6 "fast" charged tracks and then in quadrature with the 
slow pion from the D*° decay. The error on the background shape has been discussed above. 
We take a conservative estimate of the systematic error due to this source of ±11%. We 
use the current Particle Data Group values for the relevant D*° and D° branching ratios of 
(61.9±2.9)% (D*° -> n°D°) and (3.83±0.09)% (D° -> K-<k+) 0. The relative errors, 4.7% 
for the D*° branching ratio and 2.3% for the D° are added in quadrature to the background 
shape error, the 7r° detection efficiency uncertainty and the tracking error. The total positive 
systematic error is 19%. We also allow for cross-feed backgrounds amounting to 4 events 
giving a total negative systematic error of 20%. 

In Fig. |5](a) we show the D*°tc + invariant mass spectrum, obtained by fitting the number 
of B events as a function of D*°n + mass. (There are two combinations per event.) We 
also show for comparison the D*°n~ mass spectrum, where no structure is expected. We 
see evidence for an excess of events in the region of the D** + (2420) and D** + (2460). There 
are four D** mesons. Two have relatively narrow widths and decay into D*tt, whereas two 
are wide, with only one decaying into D*tt [J8J. It is difficult to quantitatively evaluate the 
fraction of total D** + production in our data. We find that ~70% of the signal has one mass 
combination between 2.3 and 2.6 GeV. 

In Fig. |6| we show the 7r + 7r~7r~ mass spectrum when the D* tt + mass is required to be 
between 2.3 and 2.6 GeV. Here we fit the B yield as a function of 7r + 7r~7r~ mass. There is 
no clear feature present. 
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FIG. 4. The invariant mass spectrum of tt + tt + tt tt for the final state D*°ir + ir + ir ir , found 
by fitting the B yield in bins of Att mass. 



Let us now see how the presence of this final state affects the LLW prediction. In Fig. |7| we 
show the CLEO data for (B — > .D* + 7r + 7r~7r~7r ) plotted as a function of the four- 
pion invariant mass squared, normalized to the semileptonic rate |7]], and compared with 
the LLW prediction 0. We also plot (B — > /}*°7r + 7r + 7r~7r^ , again normalized to the 
semileptonic rate. In principle a non-zero rate in the D*°(47r)° final state is indicative of an 
additional contribution to the D* + (4n)~ final state, beyond what is expected in factorization, 
and needs to be subtracted to make an accurate prediction.^ In fact, the D*°(47r)° rate is 
consistent with zero in the mass squared region covered by the LLW prediction. 

In conclusion we have made the first measurement of B [B — > D*°7r + 7r + 7r~7r~) = (0.30 ± 
0.07 ± 0.06)%. The reaction has a large component of D** + — > D*°tt + . We determine the 
relative rate 



Rn 



TIB 



D*° 



n n % % 



T [W 



D* + 7l + 7l~7l 7T° 



0.17 ±0.04 ±0.02 . 



(3) 



We have no evidence in the D*° final state for 4tt° masses below 2.9 GeV 2 and set an upper 
limit on i?o- < 0.13 at 90% confidence level in this restricted mass region. 

LLW have used the B — > D* + 7t + 7t~7t~tt reaction to test the 4tt mass dependence of 
factorization. They point out that a perturbative origin for factorization should cause a 



2 If the D*°(47r)° final state comes from D** + production, then some D** + , (half as much according 
to isotopic-spin symmetry) should be present in the D* + (Att)~. 
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FIG. 5. The invariant mass spectra of (a) D*°7r + , and (b) D*°7r~ for the final state 
.D* 7r + 7r + 7r~7r~, found by fitting the B yield in bins of D* ^ mass. (There are two combina- 
tions per event.) 

weakening of the prediction with increasing 4n mass. However if the basis for factorization 
is the large N c limit, where N c refers to the number of colors, no such weakening should 
occur. LLW also suggest that the presence of D** production might cause their factorization 
test to be inaccurate. We have found such a presence in the analogous reaction B — > 
.D*°7r + 7r + 7r~7r~. However, the An mass region that is populated is higher than that used by 
LLW, so no effect on their prediction can be inferred. 

We thank Z. Ligeti, M. Luke and M. Wise for interesting discussions. We gratefully 
acknowledge the effort of the CESR staff in providing us with excellent luminosity and 
running conditions. This work was supported by the National Science Foundation, the U.S. 
Department of Energy, the Research Corporation, the Natural Sciences and Engineering 
Research Council of Canada and the Texas Advanced Research Program. 
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FIG. 6. The invariant mass spectra of tt^tt it when the D*°ir + invariant mass is between 2.3 
and 2.6 GeV, found by fitting the B yield in bins of 7r + 7r~7r~ mass. 




FIG. 7. ^gj \B° — > D*TnnrTT ) , where M 2 is the 47r invariant mass squared, normalized to the 

semileptonic width T(W D*£-u). The circles are the CLEO data for S° -> Z?* vr + 7r + 7r"7r", 
the squares the model prediction of LLW and the triangles are the data for B° — > Z?* + 7r + 7r _ 7r~7r°. 
There is an additional systematic normalization uncertainty on the triangle points of 9% and on 
the circles of 16%. 
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